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We report a systematic investigation of Ba(Fei_iCo2,)2As2 based on transport and ^^As 
NMR measurements, and establish the electronic phase diagram. We demonstrate that doping 
progressively suppresses the uniform spin susceptibility and low frequency spin fluctuations. 
The optimum superconducting phase emerges at Xc — 0.08 when the tendency toward spin 
ordering completely diminishes. Our findings point toward the presence of a quantum critical 
point near Xc between the SDW (spin density wave) and superconducting phases. 
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The recent discovery of iron-pnictide high super- 
conductors^"^ poses a new intellectual challenge in con- 
densed matter physics. Although the superconducting 
mechanism remains unknown, it has become apparent 
that iron-pnictides share remarkable similarities with 
high Tc cuprates: the undoped parent phases LaFeAsO 
and BaFe2As2 have FeAs square-lattice sheets, and itin- 
erant electrons in these layers undergo antiferromagnetic 
long range order (AFLRO) below modest temperatures, 
TsDW ~ 140 K;'^"^ electron or hole doping suppresses 
the AFLRO and induces superconductvity. There exist 
clear dissimilarities, too. For example, all five 3d orbitals 
of Fe atoms contribute to the multiple Fermi surfaces, 
hence to superconductivity, in iron-pnictides. In con- 
trast, only the Cu 3d2,2_j,2 orbital plays a role in cu- 
rates. Furthermore, substitution of Co atoms into Fe sites 
of the parent phases results in electron doping and in- 
duces superconductivity^^"^^ without creating localized 
moments, while Zn^+ ions doped into Cu^+ sites in- 
duce localized moments and destroy superconductivity 
in cuprates. Sorting out these similarities and dissimilar- 
ities may lead us to an understanding the mechanism of 
high Tc superconductivity in iron-pnictides as well as in 
cuprates. 

In this Letter, we utilize transport and NMR tech- 
niques to probe the evolution of the electronic properties 
of Ba(Fei_a;Co2;)2As2 single crystals from the undoped 
SDW phase (a; = with Tsdw = 135 K), underdoped 
SDW phase {x = 0.02 with Tsdw = 100 K, and x ^ 0.04 
with Tsdw = 66 K), optimally doped superconducting 
phase {x = 0.08 with T^ = 22 K), to the shghtly over- 
doped regime {x = 0.105 with T^. = 15 K). We establish 
the electronic phase diagram which has a quantum crit- 
ical point at ~ 0.08, and explore the possible relation 
between paramagnetic spin fluctuations and the mecha- 
nism of superconductivity. From NMR Knight shift and 
spin-lattice relaxation rate measurements, we show that 
electron doping progressively suppresses the uniform spin 
susceptibility Xspin and low frequency spin fluctuations. 



The optimally doped superconducting phase emerges at 
x = 0.08 (~ Xc) as soon as doped electrons completely 
suppress the tendency toward AFLRO. 

We grew single crystals of Ba(Fei_a;Coa;)2As2 from 
FeAs flux,^^ and determined the actual Co concentra- 
tion X by EDS (Energy Dispersive X-ray Spectroscopy). 
The X = 0.08 sample is the same specimen as the one 
we previously studied and labeled as a; = 0.1,^'^' and 
has the nominal starting composition of x = 0.1. The 
bulk characterization of the x — 0.04 and 0.08 crystals is 
discussed elsewhere. We carried out NMR measure- 
ments using standard pulsed NMR techniques for aligned 
single crystals with the total mass of ~ 2 to ~ 20 mg. We 
measured the frequency dependence of the spin-lattice 
relaxation rate divided by temperature T, jYt, within a 
single NMR line by integrating the EFT signals for an ap- 
propriate frequency range. AC resistivity measurements 
were performed with the standard four probe method. 
We define the critical temperature Tsdw of AFLRO 
based on the maximum negative slope of Pab,^^ while 
Tc was determined by SQUID measurements. We sum- 
marize the resistivity data in Fig.l, and the electronic 
phase diagram of Tsdw and Tc in Fig. 2. 

In Fig. 3, we present typical ^^As NMR lineshapes of 
the Iz = —1/2 to -1-1/2 transition at 175 K obtained 
by EFT of the spin echo envelope measured in a fixed 
magetic field of Bext ^7-7 Tesla applied along the crystal 
c-axis. For clarity, we converted the horizontal axis of the 
NMR resonance frequency / to the NMR Knight shift, 
K, based on the following relation; 

f ^ {1 + KhnBc^t, (1) 

where jn is the nuclear gyromagnetic ratio, and 

^ — K spin "t" KcfiQjji^ 

Ahf 



(2a) 

Naiib^'""" ^^^^ 

The Knight shift K consists of two separate contribu- 
tions. The first term in Eq.2a, Kgpim is the spin con- 
tribution, and proportional to the local spin susceptibil- 
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ity along the c-axis, Xspin, as defined in Eq.2b. Tiie sec- 
ond term in Eq.2a, Kchem, is a temperature-independent 
background term {the chemical shift), and is not re- 
lated to Xspin- From the measurcmciiit at 4.2 K in the 
superconducting state below Tg of the x = 0.08 sam- 
ple, we estimate Kchem = 0.224 % along the c-axis. 
Afif = 18.8 kOe/fiB in Eq.2b represents the hyperfine 
interaction along the c-axis between the observed ^^As 
nuclear spin and electron spins/^ and Na is the Avo- 
gadro's number. 

We summarize the T and x dependences of the center 
of gravity of K in Fig. 4. Eq.2b implies that a change of 
the spin contribution to the Knight shift by the amount 
of AKspin = 0.0339 % corresponds to that of i^Xspin = 
1 X 10~* emu/mol-Fe along the c-axis. We show the con- 
version on the right vertical axis of Fig. 4. The magni- 
tude Xspin is comparable to that of LaFeAsOo.gFo.i.^^ 
Notice that the origin on the right axis, Xspin — 0, is 
matched with the temperature independent background 
contribution Kchem = 0.224%. Fig. 4 represents the first 
successful measurements of the intrinsic behavior of Xspin 
for a broad range of concentrations of iron-pnictides. 

We notice two striking features in Fig. 4. First, doped 
electrons suppress Xspin- We recall that doped holes en- 
hance Xspin in high Tc cuprates,^" and the enhancement 
is interpreted as a consequence of the suppression of some 
energy scale with doping, such as the effective spin-spin 
exchange interaction J or spin pseudo-gap Apc-^^'^^ If 
we interpret the results of Fig. 4 in an analogous man- 
ner, our results imply that the energy scale increases 
with doping in Ba(Fei_a;Cocc)2As2. Second, the T de- 
pendence of Xspin is qualitatively similar in the parama- 
genetic state of all samples from undoped antiferromag- 
net (x — 0) to the slightly overdoped superconductor 
{x = 0.105); 

Xspin decreases roughly linearly for all sam- 
ples below 300 K, then begins to level off below ~ 100 K 
toward Tc in the superconducting samples with x ^ 0.08. 
Pab also shows a T-linear behavior in the same temper- 
ature regime below ~ 100 K, as initially reported for 
X = 0.08.^'^''^® That is, the high Tc superconductivity of 
iron-pnictides arises from a novel electronic state with 
Xspin ~ constant and pab ~ T after Xspin is suppressed. 

What are the dominant mechanisms that govern the T 
and X dependences of Xspin^ Since the undoped sample 
undergoes AFLRO below Tsdw = 135 K, one obvious 
possibility is that the growth of antiferromagnetic short- 
range order (AFSRO) suppresses Xspin with decreasing 
T. However, AFSRO alone can't account for the x depen- 
dence, because Xspin is smaller for the Co doped samples 
without AFLRO. In other words, if we attribute T and x 
dependencies entirely to the effects of AFSRO, we would 
be led to an unphysical conclusion that the antiferromag- 
netic spin-spin correlation is stronger in the Co doped 
superconducting phase than in the undoped phase with 
AFLRO. Moreover, growth of AFSRO generally results 
in enhancement of 1/TiT with decreasing T, as observed 
for La2_a:Sr2,Cu04,^^ but our 1/TiT data contradict with 
such a scenario as explained below. While AFSRO is very 
likely to play a significant role in controlling the behavior 
of Xspin, there must be an additional mechanism which 
suppresses Xspin- We will address this issue below. 



Before proceeding, we would like to comment briefly 
on the implications of the systematic line broadening of 
the NMR lines with doping in Fig.3.^'' Eq.2b implies that 
the distribution of Kspi in Fig. 3 reflects that of Xspin- 
Accordingly, the FFT lineshapes in Fig. 3 represent a his- 
togram of the distribution of Xspin in each sample. For 
example, some parts of the x = 0.08 sample have Xspin 
as large (small) as that of the x ~ 0.04 {x = 0.105) sam- 
ple, i.e. the local electronic properties are inherently in- 
homogeneous at a microscopic level. We emphasize that 
one can't attribute the present finding to macroscopic 
inhomogeneity of the Co concentrations. As shown in 
Fig.l, the electrical resistivity pab is very sensitive to 
phase transitions at Tsdw and Tc. If the x = 0.08 crys- 
tal had a macroscopic domain with Co 4% doping, for 
example, we would observe an additional SDW anomaly 
in the pab data. We recall that analogous microscopic in- 
homogeneity was also observed for the hole-doped high 
Tc cuprates La2_a:Srj,Cu04,^^'^'* and understood as the 
consequence of microscopic patch-by-patch inhomogene- 
ity of the carrier concentration with the typical length 
scales of ^ 3 nm. Our results in Fig.3 show that anal- 
ogous "patchy" electronic inhomogeneity also exists in 
Ba(Fei_j;Coa;)2As2 with comparable length scales. 

Additional clues on the nature of spin correlations may 
be seen in the T and x dependencies of \/T\T. We sum- 
marize 1/T\T measured at the peak of the NMR line- 
shapes in Fig. 5. Theoretically, 1/TiT is related to the 
wave vector integral of the low frequency component of 
spin fluctuations, 

^~Ei7n^(q)i^^^^^' (3) 

where Im x(q, /) is the imaginary part of the dynamic 

electron spin susceptibility at the NMR frequency / 
65 MHz), |A(q)|2 = \Acos{q^a/2)cos{qya/2)\^ is the 
form factor of the transferred hyperfine interaction at the 
''^As sites (a is the lattice constant), and the wave vec- 
tor summation of q is taken over the first Brillouirn zone. 
The undoped a; = sample does not exhibit a large en- 
hancement of 1/TiT immediately above Tsdw = 135 K, 
in agreement with Kitagawa et al.^^ There are two rea- 
sons for the absence of the signature of critical slowing 
down. First, the AFLRO is accompanied by a tetrag- 
onal to orthorhombic structural phase transition,* and 
the simultaneous magnetic and structural phase transi- 
tions are weakly first order. -'^^ Second, since wc are ap- 
plying Bcxt along the c-axis, the ab-plane components 
of the dipole hyperfine fields from Fe moments, as well 
as the transferred hyperfine fields in Eq.3; cancel out at 
the location of As sites when spin correlations are com- 
mensurate'^' * with the lattice. 

On the other hand, underdoped x = 0.02 and 0.04 sam- 
ples show a divergent behavior of l/TiT due to critical 
slowing down of spin fiuctuations toward Tsdw = 100 K 
and 66 K, respectively. We also conformed that the para- 
magnetic NMR signals disappear below Tsdw- These 
results strongly suggest the SDW state, presumably in- 
commensurate with the lattice, emerges below Tsdw as 
a result of the second order phase transition. 
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We have not conducted NMR measurements in the 
concentration range x = 0.05 ~ 0.07. However, the NMR 
lineshapes of the x = 0.04 and 0.08 samples are signifi- 
cantly superposed because of the patchy inhomogeneity, 
and we have found that 1/TiT measured at the half- 
intensity position on the lower K side of the x = 0.04 
is nearly the same as l/TiT measured at the half in- 
tensity position on the higher K side of the x = 0.08 
crystals, see Fig. 3. These 1/TiT data points measured 
in the middle between the peaks oi x = 0.04 and 0.08 
are plotted in Fig. 5, labeled as the results for the effec- 
tive concentration of "6%." The enhancement of 1/TiT 
toward Tc = 22 K for "6%" and 8% samples signals the 
residual effects of low frequency AFSRO at low temper- 
atures, and we conclude that the FeAs layers are on the 
verge of forming an SDW ground state even for x < 0.08. 

In contrast, 1/TiT for the slightly overdoped x = 0.105 
levels off after monotonic decrease with T, and exhibits 
no evidence for the enhancement due to low frequency 
AFSRO. Equally important to notice is that 1/TiT ini- 
tially decreases with T below 300 K even in the under- 
doped X = 0.02 and 0.04 samples. We confirmed that 
1/TiT at the Co sites in a; = 0.08^^ and x = 0.04 
(not shown) shows qualitatively the same behavior as at 
'^^As sites; hence the present results are not the conse- 
quence of accidental cancellation of the hyperfine fields 
at ^^As sites. Recalling that the uniform (i.e. q = 0) 
spin susceptibility Xspin is also suppressed below 300 K 
for all concentrations, we conclude that the low energy 
spin excitations are suppressed with decreasing T except 
near Tsdw, i-e- Ba(Fei_xCox)2As2 exhibits spin pseudo- 
gap behavior for a broad concentration range. We note 
that analogous pseudo-gap behavior was first observed 
in LaFeAsO0.89F0.11 {Tc = 28 K) using ^^F NMR.^s Al- 
though the mechanism behind the spin pseudo-gap be- 
havior is unclear at this time, it is interesting to realize 
that the Ji — J2 Hcisenbcrg model, which is considered 
a viable starting point for the theoretical description of 
iron-pnictides,^^ is indeed expected to exhibit a gapped 
behavior in NMR properties. On the other hand, if 
we consider Fe spins as localized moments described by 
the Ji — J2 Heisenberg model, we expect that the ex- 
change narrowing effects lead to 1/Ti ^ constant.^'^ 
hence 1/TiT ~ in the intermediate temperature 

range T ~ Ji,2/2 or higher. No evidence for such lo- 
calized moment behavior is seen in Fig. 5 at least up to 
300 K. 

How does electron doping affect the spin pseudo-gap? 
In Fig. 4, we fit the temperature dependence of K by as- 
suming an empirical activation formula, K = A -\- B x 
exp{—ApG/kBT), where A and B are constants. The 
best fit yields Apa/ke 710 K, 570 K, 520 K, and 490 K 
for X = 0, 0.04, 0.08, and 0.105. It is not clear whether 
the apparent systematic decrease of Apc with doping is 
intrinsic, because we found that the fit is extremely sen- 
sitive to the choice of the constant term A and the uncer- 
tainty of ApG is as large as ~ 30 %. For LaFeAsOi-xF^;, 
analogous fit of 1/TiT or yields Apc/kB ^ 160 K for 
x>0.10.9'i9>28 

In order to gain insights into the essential physics of 
the suppression of low energy spin excitations without 



making any assumptions, we create the contour map of 
equal values of K or 1/TiT in Fig. 2. First, let us focus 
our attention on the black dotted lines representing the 
equal value contours of K. We can clearly observe, with- 
out relying on any theoretical models or assumptions, 
that Xspin is suppressed for lower T and larger x. Turning 
our attention to the blue dashed lines representing the 
equal value contours of 1/TiT, we also observe that the 
residual tendency toward forming the SDW ground state 
barely disappears at the optimum doping x = 0.08 with 
the highest Tc = 22 K. Thus we conclude that the "su- 
perconducting dome" appears right on top of the quan- 
tum critical point Xc ~ 0.08 where the tendency toward 
forming the SDW phase disappears. We note that the 
disappearance of the low frequency components of spin 
fluctuations does not mean that spin fluctuations dimin- 
ish all together. In fact, the pseudo-gap behavior suggests 
that the over all spectral weight of spin fluctuations may 
be shifted to higher energies. 
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Figure Captions 
Fig.l 

The in- plane resistivity pab- Vertical solid arrows mark 



Fig.2 

The electronic phase diagram of Ba(Fei_a;COa;)2As2. 
The filled circles and triangles represent Tsdw and Tc, 
respectively. The blue dashed lines and black dotted 



lines represent the equal magnitude contours of 1/T\T 
and the Knight shift K, respectively. 

Fig.3 

Representative ^^As NMR lineshapes obtained by FFT 
at 175 K. We converted the horizontal frequency axis to 
that of K using Eq.2. Also plotted is l/TiT measured 
at various points within each FFT line. The horizontal 
bars specify the range of K in which the FFT intensity 
was integrated for 1 /T\T measurements. 

Fig.4 

The center of gravity of the '^^As NMR Knight shift 
K measured along the crystal c-axis. K = 0.224 % 
corresponds to the background arising from the chemical 
shift Kchem- The conversion to spin susceptibility Xspin 
based on Eq.2b is shown on the right vertical axis. 
Vertical solid and dashed arrows mark Tsdw and Tc, 
respectively. Solid curves are a fit to an activation for- 
mula with Apc/kB = 711 K (a; = 0), 570 K (a; = 0.04), 
520 K{x = 0.08), and 490 K (a; = 0.105). 

Fig.5 

Filled symbols show 1/TiT measured at the peak of the 
As NMR lineshapes for 0, 2, 4, 8, and 10.5% doped 
samples. Orange open circles (triangles) labeled as " 6%" 
were deduced from the half intensity position of the FFT 
lineshape of the 8% (4%) doped sample on the higher 
(lower) frequency side. See the main text for details. 



